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Abstract 

At low temperatures and 3 He concentrations below ~ 6.6%, there is ex- 
perimental evidence about the existence in liquid helium mixtures, of stable 
vortices with 3 He-rich cores. When the system is either supersaturated or 
submitted to a tensile strength, vortices lose stability becoming metastable 
and eventually completely unstable, so that their cores freely expand. Within 
a density functional approach, we have determined the pressure- 3 He concen- 
tration curve along which this instability appears at zero temperature. 
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The structure, stability and dynamics of quantized vortices in superfluid 4 He have been 
extensively studied either experimental or theoretically (for a systematic review see and 
for recent work, see for example |2|-|6| and Refs. therein). However, lesser work has been 
done in the case of 3 He- 4 He solutions An interesting aspect of the vortex structure 
in these mixtures was disclosed by Williams and Packard J7J, who provided experimental 
evidence that, in diluted 3 He- 4 He solutions at low temperature (T), stable vortices present 
3 He condensation onto the core. 

In this work we address the problem of the stability of a vortex line as a function of 
pressure (P) and 3 He concentration (x). In particular, we will discuss the implications that 
the existence of 3 He-rich vortices may have on the critical supersaturation of isotopic helium 
mixtures, and will determine the T=0 vortex spinodal line as a function of x. 

The interest in studying the stability of vortex lines at negative pressures stems from 
recent attempts to describe the phenomenon of quantum cavitation in liquid 4 He 0. In 
particular, it has been found @ that at T=0, vortex lines become unstable at P ~ -8 bar, 
whereas the spinodal point is at P~ -9.4 bar ||, thus quantitatively showing that vorticity 
raises the spinodal line (see also ||). 

In 3 He- 4 He mixtures the situation is more complex. Indeed, a vortex may become unsta- 
ble either increasing the 3 He concentration, or submitting the solution to a tensile strength 
that may originate a negative pressure. A characteristic of the zero temperature P-x phase 
diagram that makes more intrincate the study of vortex lines in this system, is the existence 
of a demixing line Pd(x), experimentally determined from saturation up to 20 bar |J. A 
recent calculation [|10| has found that line to continue down to x ~ 2.4 % and P~ -3.1 
bar, which is the T=0 spinodal point of pure 3 He. The existence in the negative pressure, 
metastable region, of this equilibrium line between pure 3 He and the mixture, affects the 



description of cavitation caused either by bubble growing [10,11] or by vortex destabilization. 



Instabilities caused by supersaturation were studied in [F^] within the so-called Hollow 
Core Model (HCM), adapted to helium mixtures replacing the hollow by a 3 He-rich core. 
This study was carried out at positive pressures, where the model works better. We have 



predicted that at P=0, vortices become completely unstable for x >8.2%. 

Although HCM might seem too crude a model, at positive pressures and close to the 
demixing line it bears the basic physical ingredients, so it can be used as an useful guide 
to understand the appearance of unstable vortices. Let us just recall that within HCM, the 
total energy per unit vortex length as a function of the core radius R reads ||12||: 



n HCM = SR + VR 2 + E ln(Roo/R) ■ (1) 

It is the sum of an S-" surface term", plus a V-" volume term", plus a kinetic energy term. 
When the vortex is in the stable phase, all three terms are positive, and only stable vortices 
may exist. 

On the contrary, in the metastable phase the factor V in the volume term becomes 
negative. Depending on whether the system is underpressured or supersaturated, this factor 
is proportional either to the difference AP between the pressure in the hollow and in the 
bulk, or to the difference A/x between the chemical potential of 3 He in the mixture and of 
pure 3 He [IT]. As a consequence, the stable vortex becomes metastable and there exists a 



critical vortex configuration for which the potential barrier has a maximum located at a core 
radius R c . At the saddle point, when the energy difference between critical and metastable 
vortices vanishes, the vortex core freely expands. 

These arguments can be made quantitative within the density functional approach. To 
this end, we resort to the density functional of Refs. [|K],[ll]] which describes the basic ther- 
modynamical properties of liquid helium mixtures at zero temperature. For the sake of 
simplicity, we address the problem of a vortex line. Using cylindrical coordinates and taking 
the vortex line as z-axis, the density profiles depend on the r-distance to the z-axis and are 
obtained solving the Euler-Lagrange equations 

^^ = 0,, = 3,4, (2) 

where cu(p3, p±) is the grand potential density functional flT0| , |TT] | to which we have added 
a centrifugal term h 2 p^n 2 / ftm^r 2 ) associated with the superfluid flow. We choose the 



quantum circulation number n=l because it corresponds to the most stable vortex 
is the 4 He atomic mass, and p q are the particle densities of each helium isotope. 

For given P and x, Eqs. (0) are solved imposing that at long distances from the z-axis, 
p q equals that of the metastable, homogeneous liquid p q (x is simply p^/(p^ + P4)), and that 
P4 and the r-derivative of p% are zero on the z-axis. Notice that the metastable and critical 
configurations are solutions of these equations for the same P and x conditions. The barrier 
height per unit vortex length is: 

AQ = 2nJ rdr [u(p c 3 , pi) - u(jg, p™)] , (3) 

where p c q and p™ are the particle densities of the critical and metastable vortices, respectively. 
It is worth to note that AQ is a finite quantity: there is no need to introduce any r-cutoff 
as it would have been unavoidable if we had described either configuration separately. 

Depending on the situation of the metastable vortex in the P-x plane, there may exist 
two different kinds of critical configurations. To illustrate it, we show in Fig. 1 the P-x phase 
diagram at T=0 [|Tt]] and three selected metastable configurations labeled 1 to 3. The grey 
zone represents the stable region, and the dashed line is the demixing line. Configuration 1 
is underpressed, configuration 2 is supersaturated and configuration 3 is both. In all three 
cases, the metastable configuration corresponds to a rather compact vortex filled with 3 He 
whose radius increases with increasing x. This is not the case for the critical vortex. Indeed, 
as configuration 2 is in the supersaturated region, the critical vortex may have a large R c 
radius if that point is close enough to the demixing line, and its core is filled with almost pure 
3 He. R c diverges at the demixing line, and also AQ. Since configuration 1 is underpressed 
and undersaturated, the critical vortex also has a large radius provided point 1 is close 
to the P=0 line, but its core is almost empty, with the surface covered by 3 He (Andreev 
states). R c diverges at the P=0 line, and also AQ. As configuration 3 is underpressed and 
supersaturated, it has two possible critical configurations, one bearing the characteristics of 
configuration 1, and another bearing those of configuration 2. 

These possibilities are displayed in Fig. 2. Fig. 2 (a) corresponds to a type 1 configu- 



ration with P=-1.66 bar, x =1%, whereas Fig. 2 (b) corresponds to a type 2 configuration 
with P=0.91 bar, x =8%. 

Fig. 3 shows the barrier height per unit length as a function of P for x =1 to 9%. 
For the sake of illustration, we display for x =4%, the barriers corresponding to the two 
kinds of critical vortices already discussed: the dashed (solid) line is AQ for empty- (filled-) 
core configurations. Notice that these curves have different slopes because they diverge at 
different pressures, the former at P=0, and the later at P=Pd- For a given x, the P-value 
at which AQ is negligible defines a point along the vortex spinodal curve. That curve is the 
solid line in Fig. 1. 

Fig. 4 shows the core radius of the saddle configurations as a function of x (solid line). 
Following ||, we have defined that radius as the r- value at which the superfluid circulation 
current p%(r)/i has a maximum. We have found that metastable vortices in the mixture 
have a core radius larger than in pure 4 He 0. This is in agreement with the experimental 
findings for stable vortices Also shown in that figure is the radius of the stable vortex 
at P=0 (dashed line), which is actually metastable above x =6.6%. 

The above results have implications on the critical supersaturation degree Ax cr of iso- 
topic helium solutions at low temperatures. Recent experiments [ ]T^ , [Tr5] ] have found Ax cr 



below ~1%, whereas classical nucleation theory yields ~10% [0,0- The microscopically 



calculated spinodal line [17j is about the same value. Within the HCM, we have argued [12| 



that the rather small degree of supersaturation experimentally found could be due to desta- 
bilization of 3 He-rich vortices. The density functional approach yields Ax cr values around 
2% (we recall that the maximum solubility at P=0 is ~ 6.6 % |18j). This is shown in Fig. 
5 for P=0, 0.5 and 1 bar. A discrepancy with experiment still exists. It is unclear whether 
considering more realistic vortex geometries, like vortex rings, could bring theory closer to 
experiment. Other possibilities to improve on the agreement, such as vortex destabiliza- 
tion due to quantum tunneling through, or to thermal activation over the barrier, seem to 
be ruled out. The former because of the extremelly small quantum-to-thermal crossover 



temperature ||19| , and the later because of the large mass of the critical vortex [12 



Vortex destabilization at negative pressures may also have implications on the phe- 
nomenon of cavitation in isotopic helium solutions JlQ| , |TT|j . Realistic calculations |2(J in- 
dicate that for a 3 He concentration as small as x cavitation driven either by quantum 
or by thermal fluctuations triggers phase separation if the system is submitted to a tensile 
strength of 8.2 bar. If there are metastable vortices in the mixture, the present calculations 
show that for much smaller tensile strengths, of about 5.2 bar, the solution undergoes phase 
separation. A delicate question is whether 3 He atoms have enough time to diffuse into the 
vortex core on the time scale of current cavitation experiments [^TJ when their concentration 
is too small. 

This work has been supported in part by the CICYT, by the Generalitat de Catalunya 
Visiting Professors and GRQ94-1022 programs, and by the CONICET (Argentine) Grant 
No. PID 97/93. 
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FIGURES 

FIG. 1. P-x phase diagram at T=0. The grey zone represents the stable region. The dashed 
curve is the demixing line, which ends at the P=-3.12 bar, x = 2.43% cross, and the solid curve is 
the vortex spinodal line. 

FIG. 2. Panel (a): vortex profiles for x = 1% and P=-1.66 bar. Panel (b): vortex profiles for 
x = 8% and P=0.91 bar. The solid lines represent the total particle density, and the dash-dotted 
(dashed) lines, the p± (pz) densities. Critical (metastable) configurations are denoted as p c (p m ). 

FIG. 3. Barrier height per unit vortex length as a function of P for the indicated 3 He concen- 
trations. At x = 4%, the dashed line corresponds to empty-core configurations, whereas the solid 
line corresponds to filled-core ones. 

FIG. 4. Solid line, radius of the saddle vortex core. Dashed line, radius of the stable vortex 
at P=0. 

FIG. 5. Barrier height per unit vortex length as a function of x for P=0, 0.5 and 1 bar. The 
corresponding critical x-values are 8.46, 8.74 and 8.97%, respectively. 
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